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In this paper we study the nilpotent 2-step, corank 2 sub-Riemannian metrics that are 

' nilpotent approximations of general sub-Riemannian metrics. We exhibit optimal syntheses for 

^> . these problems. It turns out that in general the cut time is not equal to the first conjugate 

. ' time but has a simple explicit expression. As a byproduct of this study we get some smoothness 

1^ . properties of the spherical Hausdorff measure in the case of a generic 6 dimensional, 2-stcp 

' corank 2 sub-Riemannian metric. 



X 



Abstract 



1 Introduction 

(N 

> . 

\^ , 1.1 Sub-Riemannian manifolds 

' In this paper, by a sub-Riemannian manifold we mean a triple (M, A, g), where M is a connected 

lO ■ orientable smooth manifold of dimension n, A is a smooth vector distribution of constant rank 

lO . m < n, satisfying the Hormander condition and g is an Euclidean structure over A. 

\ A sub-Riemannian manifold has the natural structure of a metric space, where the distance is 

the so called Carnot-Caratheodory distance 



d{qo, Qi) = ^^i{j ^JS'y{t){iit),ii't)) dt \ : [0,T] ^ M IS a. Lipschitz curve, 

7(0) = qoniT) = qu W) ^ A^(t) a.e. in [0,r]}. 

As a consequence of the Hormander condition d is actually a distance inducing the topology of the 



manifold. This is the Rashevsky-Chow Theorem, see for instance [R| or [18|. 



Along this paper we assume that the structure is 2-step bracket generating i.e. 

TgM = Ag + [A, A]g, foT cvcry q £ M, 

and we quote a 2-step sub-Riemannian metric by its rank and its dimension, i.e. with the pair 
(m, n). The quantity k = n — m\s called the corank of the structure. 



This research has been supported by the European Research Council, ERG StG 2009 "GeCoMethods" , contract 
number 239748, by the ANR Project GCM, program "Blanche", project number NT09-504490 and by the DIGITEO 
project CONGEO. 
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It is know from Mitchell |2C] that the Hausdorff dimension of M, as a metric space, is f) = 
2n — m > n. In this paper we focus on the case (m, m + 2). 

A sub-Riemannian manifold is left-invariant if M = G, a Lie group, and both A and g are 
left-invariant over G. 

Locally, the pair (A, g) can be specified by the data of a set of m smooth vector fields spanning 
A, being an orthonormal frame for g, i.e. 

Aq = span{Xi{q),...,Xm{q)}, Sq{Xi{q),Xj{q)) = 6ij. (1) 

In this case, the set {Xi, . . . , X^} is called a local orthonormal frame for the sub-Riemannian 
metric. 

The sub-Riemannian metric can also be expressed locally in "control form" as follows. We 
consider the control system, 

m 

q = Y,UiXi{q), UiGM, (2) 

i=l 

and the problem of finding the shortest curve minimizing that joins two fixed points q^, (71 G M is 
naturally formulated as the optimal control problem, 



\ 



uKt) dt ^ min, q{0) = qq, q(T) = q^. (3) 



i=l 



A geodesic for the sub-Riemannian metric is a curve, parametrized by constant velocity, such that 
every short enough piece of it is a local minimizer of the length. For the sub-Riemannian metrics 
given in control form all the geodesies can be computed with Pontryagin's maximum principle |22|. 
In the 2-step bracket generating case, it is known that there is no strict abnormal minimizer, and 
all geodesies are projections on M of the trajectories of the Hamiltonian system associated with 
the following Hamiltonian over T*M 



^ m 

H{X,q) = -Y,{^,X,{q)f, XeT*M. (4) 
^ i=i 

and corresponding to the level set {H = c}, for c > 0. 



1.2 Nilpotent approximation 

Consider a sub-Riemannian manifold (M, A, g) and fix a point q £ M. The Lie bracket induces a 
skew symmetric tensor bilinear mapping 

[•,•], : A, X A, ^ r,M/A,. (5) 

Then, for every Z* G (T^M/Ag)*, we have 

Z*{[X,Y] + Aq) = {Az^iX),Y)^, 

for some g-skew symmetric endomorphism Az* of A^. 

Remark 1 (Notation) . We denote by Cq the /c-dimensional space of skew symmetric endomorphisms 
of Aq obtained by taking the union of all the Az* at q. This notation is used in the Appendix. 
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The space Lg = © TqM / Ag is endowed with the structure of a 2-step nilpotent Lie-algebra 
by setting 

[{VuWi),{V2,W2)\ = (o,[yi,y2] + Ag). 

The associated simply connected nilpotent Lie group is denoted by Gq and the exponential mapping 
Exp: Lq — > Gq is one to one and onto. By translation, the metric gg over Ag allows to define a 
left-invariant sub-Riemannian metric over Gq. 

Definition 2. The sub-Riemannian metric on Gq defined above is called the nilpotent approxima- 
tion of (Af , A, g) at q. 

Any k dimensional vector sub-space Vg of TqM, transversal to Ag allows to identify Lq and Gq 

to TqM ~ Aq®TqM/Aq. 

Fix qq G M. We can chose coordinates x in Agg such that the metric ggg is the standard 
Euclidean metric, and for any linear coordinate system y in Vgo, there are skew symmetric matrices 
Li, . . . ,Lk € so(m) such that the mapping iM) writes 



[X,Y] + A 



/ X'LiY 



\ X'LkY 



where X' denotes the transpose of the vector X. Then the nilpotent approximation written in 
control form is 



Xi — Ui, i — 1, . . . , 7/1, 

yi = ^x'Liu, 



(6) 



= \x'LkU. 

The construction of the nilpotent approximation given in Definition |2| makes sense for any 
sub-Riemannian metric, but it coincides with the standard one (see 0, ^) in the 2-step bracket 
generating case only. 

Proposition 3. The distribution is 2-step bracket generating if and only if the endomorphisms 
o/ Ag, Li, i = l,...,k (respectively the matrices Li when coordinates y in Vg are chosen) are 
independent. 

In the 2-step bracket generating case these linear coordinates y in TqM/Aq may be chosen in 
such a way that the endomorphisms Li, i = 1, . . . ,k are orthonormal with respect to the Hilbert- 
Schmidt norm {Li,Lj) = ^trace{L^Lj). This choice defines a canonical Euclidean structure in 
TqM/Aq and a corresponding volume in TqM/Aq. Then using the Euclidean structure over Ag 
we get a canonical Euclidean structure over Ag © TqM/Aq. The choice of the vector subspace 
Vg induces an Euclidean structure on TqM which depends on the choice of Vg, but the associated 
volume on TqM is independent on this choice. 

Definition 4. This volume form on M is called the Popp measure. 
The Popp measure is a smooth volume form. 
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1.3 Statement of the results 



1.3.1 History 

The main purpose of this paper is to build the optimal synthesis for (m, m + 2) nilpotent sub- 
Riemannian metrics, i.e. the set of all trajectories starting from the identity of the group and 
realizing the minimum of the distance, with a precise description of their cut time. 

Optimal syntheses are in general very difficult to obtain. Usually the steps are the following: 

- Apply first order necessary conditions for optimality (which in the case of sub-Riemannian 
manifolds are given by the Pontryagin Maximum Principle) to reduce the set of candidate 
optimal trajectories. This first step can be already very difficult since one should find solutions 
of a Hamiltonian system, which is not integrable in general. 

- Use higher order necessary conditions to reduce further the set of optimal trajectories. This 
step usually leads to the computation of the conjugate locus, i.e. the set of points up to which 
geodesies are locally optimal. 

- Prove that no strict abnormal extremal is optimal (for instance by using conditions such as 
the so called Goh condition 0, §]). If one fails to go beyond this step, then one can hardly 
get an optimal synthesis, since no general technique exists to treat abnormal minimizers. 

- Among all solutions of the first order necessary conditions, find the optimal ones. One has to 
prove that, for each point of a candidate optimal trajectory, there is no other trajectory among 
the selected ones, reaching that point. The first point after which a first order trajectory loses 
global optimality is called a cut point. The union of all cut points is the cut locus. 

As a consequence of these difficulties, optimal syntheses in sub-Riemannian geometry have been 
obtained in few cases. 

The most studied cases are those of left invariant sub-Riemannian metrics (see for instance Q 
for a classification in the 3D case). The first optimal synthesis was obtained for the Heisenberg 
group in ||l6| , [l^. Then complete optimal syntheses were obtained for the 3D simple Lie groups 
SU{2), SO{3), SL{2), with the metric induced by the Killing form in 12 1. An impressive work 



has been done by Yuri Sachkov who obtained the optimal synthesis for the group of motions of the 
plane 5^(2) (see ||). 

In dimension larger than 3, only nilpotent groups have been attacked. The complete optimal 
synthesis was obtained in H in the contact nilpotent case. Some results were obtained by Y. 



Sachkov for the Engel and Cartan groups [2^, |24 1 . 



When a Lie group structure is not available there are also some results: the optimal synthesis 
was obtained for a neighborhood of the starting point in the 3D contact case in [^, 0, 15 1 and in 



the 4D quasi-contact case in |14]. The optimal synthesis was obtained in the important Martinet 
nilpotent case, where abnormal minimizers can be optimal (see ||5|). They also solved the problem 
for certain perturbations of this case where strictly abnormal minimizers occur (see p!o[| ). 
To our knowledge, no other case has been solved. 

It is interesting to notice that when the sub-Riemannian metric is invariant by certain continuous 
transformation (e.g. rotations) then most of the cut points are automatically conjugate. This 
happens for instance on the Heisenberg group, on SU (2) and in the contact nilpotent case. 
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Remark 5 (Notation). In the case of our nilpotent approximations, covectors in T*M can be iden- 
tified with vectors in TqM via the Euchdean structure of TqM given by the choice of Vq. In our 
coordinates these covectors/ vectors are typicahy denoted by {uo,r). 

For nilpotent (m, m + 1) sub-Riemannian metrics that are nilpotent approximations of general 
sub-Riemannian metrics, the control systems can be written as 



Ui, i = l,...,m, 
^x'Lu, L skew s} 

Denote by ct(L) the set of all moduli of eigenvalues of the matrix L. The following fact is proved 



i = i.'Lu. L skew symmetric, 



m 



Theorem 6. (j^) Arclength geodesies of system ^Cl\} , starting from the origin, are parametrized 
by an initial covector Aq = {uo,r) € S"^~^ x W, and they are optimal until time 



2tt 





= Ui, i 








= \x'L2U, 



\r \ maxcj(L) ' 

with the understanding tcuti^o) = +oo ifr = 0. Moreover tcut(Ao) = tconj(Ao)- 

The proof of this result is based on the fact that geodesies can be expressed in terms of usual 
trigonometric functions and, thanks to a certain monotonicity property, the cut locus can be ex- 
plicitly computed and is exactly equal to the conjugate locus. 

1.3.2 Optimal synthesis for the nilpotent {m,m + 2) case 

The main result of this paper is the optimal synthesis in the case of a nilpotent approximation in the 
{m,m + 2) case. In this case the control system can be written in coordinates q = (xi, . . . , Xk,yi,y2) 
as 

1,... ,m, 

(C2) 

Li,L2 skew symmetric. 
Set ri = |r|cos0, r2 = \r\sm6, and Lg = cos{6)Li + sm{6)L2. 

Theorem 7. Arclength geodesies of system (|C2|), starting from the origin, are parametrized by an 
initial covector Aq = {uo,r) S x M?, and they are optimal until time 

2tt _ 2-K 

tcut[M) — 7 — -p — ; — — -j— j rrT' 

inaxa[riLi + r2L2) |r| max(T(-L0j 

with the understanding tcuti^o) = +oo if r = 0. Moreover, in general, tcuti^o) 7^ iconj(Ao)- 

The last statement in the Theorem ^ says that in the corank 2 case the cut and the conjugate 
time coincide only in some particular cases, which we explicitly describe in the (4, 6) case (see 
Theorem ^) . 

The reason why the corank 2 case is more difficult than the corank 1 case is precisely the fact 
that the cut locus is not equal to the conjugate locus. (The latter we are not able to compute 
explicitly.) 

Explicit expression of geodesies for this optimal synthesis are given in Section M. 
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1.3.3 The nilpotent (4,6) case 

In the nilpotent (4, 6) case our first result is the following: 
Theorem 8. The following properties are equivalent: 
(PI) The first conjugate locus is equal to the cut locus. 

(P2) The linear coordinates y in TgM/Ag can be chosen in such a way that the pair {Li,L2) of 
4x4 skew symmetric matrices belongs to the set {Q L) Q)'^ . 

Here Q (resp. Q) denotes the set of pure quaternions (resp. pure skew quaternions), see 



Appendix 4.1 



Our second result is a continuation of the paper Q for corank 1, where the following result is 
proved for general sub-Riemannian metrics. 

Theorem 9 ([Q]). In the {m,m + l) case the Radon-Nykodym derivative of the spherical Hausdorff 
measure with respect to the Popp measure is a function, but is not in general. 

Here we show the following result 

Theorem 10. For a generic (4, 6) sub-Riemannian metri(^, the Radon-Nykodym derivative of the 
spherical Hausdorff measure with respect to the Popp measure is . 

In the previous paper Q it is shown that the Radon-Nikodym derivative of the spherical Haus- 
dorff measure with respect to the Popp measure is inversely proportional (as a function of q) to the 
volume of the unit sub-Riemannian ball of the nilpotent approximation at q. Then Theorem |l^ is 
a byproduct of the optimal synthesis given here. 

Note that in the corank 1 case, the higher differentiability of the Radon-Nikodym derivative is 
due to the fact that the conjugate locus is equal to the cut locus, which is not the case here. 

Due to the complexity of the computations even in this low dimensional case, it is not easy to 
determine the real degree of differentiability of Hausdorff measure. This is still an interesting open 
question. 

1.4 Organization of the paper 

Section |2| is devoted to the construction of the optimal synthesis for (?n, m -\- 2) nilpotent sub- 
Riemannian metrics and, as a consequence, to the proof of Theorem ^. In Sections |2.1| and 2.2 



we 



compute the exponential map. In Section 2.3 we prove that geodesies are optimal up to tcut- Finally 
in Section 2.4 we show that the cut time does not coincide, in general, with the first conjugate time. 
In Section ^ we give the proofs of Theorems ^ and 

In the Appendix we recall basic facts about quaternions, we prove a technical Lemma, and 
applying an Abraham's transversality theorem, we prove that, generically, for the (4, 6) case, a 
certain "bad set" is made of isolated points, which permits to conclude about the differentiability 
of the Radon-Nikodym derivative (Theorem p^). 



^which means for an open and dense subset of all (4, 6) sub-Riemannian metrics, endowed with the Whitney 
topology. 
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2 Exponential map and synthesis 
2.1 Hamiltonian equations in the {m,n) case 

The purpose of this section is to compute the exponential map, i.e. the set of all geodesies, 
parametrized by length, starting form the origin of the control system (^, i.e. the system 

ixi=Ui, i = l,...,m, 

\yh = ^x'LhU, h = l,...,k. 

Let Lfi = (b'lj), for h = l,...,k. Then the control system can be written in the form q = 
YlT^i UiXi{q) where q = {x,y) and 

Xi = d^^ + ]^Y^ b'ljXjdy^ , i = l,...,m. 
Setting Yfi = dy^, for h = 1, . . . ,k, the commutation relations are 

k 

[Xi , Xj] =Y,h'i^Yh, i,j = l,...,m, (8) 

h=l 

[Xi,Yj] = [Yj,Yh]=0, i = l,...,7n, j,h = l,...,k. (9) 

Define the functions on T*M, that are linear on fibers, 

Ui{X, q) = {X, Xi{q)) , z = l,...,m, (10) 
rh{X,q) = {\,Yh{q)), h = l,...,k. (11) 

These functions can be treated as coordinates on the fibers of T*M to solve the Hamiltonian system 
given by the Pontryagin Maximum Principle, see Section LI. This Hamiltonian system is associated 
with the Hamiltonian 

^ m ^ m 

Hi>^, 9) = 2 E(^' ^^(«))' = 2 E ^ ^ (12) 

i=l i=l 

Remark 11. The geodesies parametrized by length correspond to the level set {H = 1/2}. Notice 
that, for systems of type q = YlT^i UiXi{q), with fixed initial and final points, the problem of finding 
length-parametrized curves minimizing the length, is equivalent to the problem of minimizing time 
with the constraint {||m|| < 1}. 

For a function a € C°°{T* M) we have that, along the sub-Riemannian flow 

m 

d = {a, H} = ^{a, ujuj, (13) 

i=l 

where {a, b} denotes the Poisson bracket of two functions in T*M. The following Lemma gives a 
way of computing the covector A(t), solution of the Hamiltonian system associated with (|l^) in the 
coordinates {u, r). 
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Lemma 12. If u{t) and r{t) are solution of (13) corresponding to level set {H = 1/2}, then they 
satisfy 

(u{t) = {riLi + . . . + rkLk)u{t), u{0) = uq, ||uo|| = 1, 
[r{t) = 0. 

Proof. Remind that, if aj(A, q) = (A, Zi{q)), for some vector fields Zi, i = 1,2, then 

{01,02} = (A, [Zi,Z2]) . 
Applying (|^) for a = r^ and using (P) we get 

A: 

fh = ^^{rh, Ui}ui =0 rh = const. 



h=i 



Similarly, using (11), one find 



^{ui,Uj}uj = ^h^^.rki 



1=1 



□ 



Remark 13. In the following geodesies are parametrized by the initial covector A(0) = (p(0), r(0)) 
(uq, r), since r = const and Ui = (A, Xj) and at the starting point we have Xj(0) = dxi- 



2.2 Exponential map in the corank 2 case 

From now on we focus on the case {m,m + 2), i.e. when the corank k is equal to 2. We can write 
the equation of geodesies starting from the origin as follows 



s(nLi+r2L2)^^^g^ X(0) = 0, 



x{t) 



yiit) = 77 / <syLi uis)ds, yi(0) = 0, 
y2{t) = 1 [ x{syL2u{s)ds, 2/2(0) = 0. 



(14) 



Remark 14 (Notation). In the following we denote by E'^^^''^^"^ [t) the geodesic, parametrized by 
the length, and starting from the origin, defined by equations ([T^), associated with Li,L2. 

Definition 15. The matrices Li,L2 being fixed, the exponential map is the map £ : x A — ?> M" 
defined by 

£{t,uo,ri,r2) = ^ = {K,ri,r2),no € 5™-\r, G E} 



Remark 16. The optimal control problem 



Xi 


= Ui, 


i 


' yi 


= \x'Liu, 






= \x'L2U, 





i = 1, . . . ,m, 



(15) 



is invariant with respect to the following change of coordinates 
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(a) orthogonal changes of coordinates in the x space, 

(b) hnear changes of coordinates in the y space. 

Indeed, let M be a nonsingular orthogonal matrix {M~^ = M') and define the new coordinates 
X = Mx. Then 

X = Mx = Mu =: u, 

and 

Vi = x'LiU = {Mx)'MLiM'{Mu) = x'MUM'u. 

Hence, in the new coordinates, Li is changed for Li := MLiM'. 
Also, it is easy to see that the change of coordinates 

yi = aiyi + 02^2, m = /3iyi + /32?/2, 

corresponds to the change 

Li = aiLi + a2iv2, L2 = fiiLi + [52L2. 

In other words we can change Li and L2 up to congruence and linear combinations. 
Using these arguments one immediately gets 



(16) 



Lemma 17. Let (ri, r2) =: (r cos 0, r smO) and Lg := cos OLi + sm6L2, Lg := — sinOLi + cos 0L2. 
Consider the rotation matrix Rg = -sirfe cose) '^'^^ orthogonal matrix M such that MLgM' is 
block diagonal. 

Denote = [^q jl^) and uq = Muq. We have the equality 



(17) 



Thanks to Lemma 17, one can always restrict to geodesies of the type 

rt 

e'''^<>uods, 
-t 



x{t) -- 

1 

yi{t) = o x{s)'Lgu{s)ds, 
^ Jo. 



y2{t) 



where Lg is in the block-diagonal form 



/ ai 

-ai 



Lg 



\ 



a 



x{s)' Lg u{s)ds, 



or 



/ ai 

-ai 



v 



ae 
-ae 



0/ 



(18) 



depending on the fact that m is even (m = 2i) or odd (m = 2^ + 1), and where the geodesic is 
associated with the covector (ri,r2) = (r, 0). 
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Remark 18. When we deal with a fixed sub-Riemannian metric we can assume also that the coor- 
dinates in the x space are chosen in such a way that ai > Oj, for every i. In this case 

2tt 27r 



rmax((j(Le)) air 
2.3 Computation of the cut time 

In this section we prove Theorem]^, i.e. we compute the last time at which a geodesic parametrized 
by length is optimal. 

We first consider the case r = 0. In this case equations ( |l4|) can be easily integrated and gives 
the straight lines 

\x{t) = Uot, 

\y^{t) = o. 

This trajectory is optimal for any time (i.e. tcut = +oo) since the sub-Riemannian length of a 
geodesic coincides with the Euclidean length of its projection on the horizontal subspace (xi, . . . , Xm), 
as follows from formula (|3|). 

In what follows we use the notation A = rLg, A = rLg and we focus on the case when A is 
even dimensional (i.e. m = 2i) and invertible. The case A non invertible (in particular A odd 
dimensional) needs an obvious modification of the proof. 



With this notation the system (18) is rewritten as 



x{t) = A-\e^^ - I)uo, 

2/1 W = 2 / x{syAu{s)ds, (19) 

If* ~ 
y2(i) = - / x{syAu{s)ds, 



2 JO 

2.3.1 Maxvi^ell points 

Consider 7(t) = £'"'''^''^(t), the geodesic associated with the problem (p^) and with initial covector 
{uq, r, 0), r > 0. Let us first show that there exists another geodesic reaching the point 7(1"*) in time 
r* = 27r/(air). Using Arnol'd's terminology, points reached in the same time by more than one 
geodesic are called Maxwell points. At the end of this section we prove that 7 cannot be optimal 
after T*. 

Set uq = (ui, n2, Ua, . . . , Um) and consider the following variation of the horizontal covector 
Uq = (cosw ui + sincj U2, — sincj ui -|- coso; it2, na, . . . , Um), ^ G [0, 2tt]. 
Denote 7'^(i) = {x'^ ,yi ,1/2) '■= E^o'^'^it) the geodesic associated with this variation. 
Claim: There exists a; / such that 7(r*) = 7^(r*). 
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Proof of the Claim. Denote by M^it) := A ^(e*"^ — I) and notice that 



MA{t) = A-\e'^ - I) 



In other words we can write 



sin airt 1— cos airt 
air air 

~l+cos airt sin airt 
air air 



\ 



sin aert 1— cos aert 

air air 

-1+cos airt sin airt 

air air 



MA{t) 



my 



where 



Di{t) 



sin ajrt 1— cos g^rt 



air 



air 



-1+cos ajrt sin ajrf 



air 



^s\ii{airt/2) / cos(ajrt/2) sin(ajrt/2)\ 
sin(ajrt/2) cos{airt/2) J ' 



a,;r 



We prove our claim by steps. 

(i). From (pO|) it is easy to see that 



smce e 



x^{T*) - x{T*) = A-\e^''^ - I)iu^ - uq) = Ma{T*){u^^ - uq) = 0, Va; € [0,27r], 
'^*^ — I (and so Ma{T*)) has its first 2x2 block equal to zero. 



(ii). Now we show that yf{T*) = yi(T*) for all uj G [0,27r]. Indeed from we get 



yi{t) 



-sA 



1 /■* 

-u'q / (e^^-/)fisno 
^ Jo 

1 1 

- (M^(t)uo,no) + -t||uo|p, 



and 



(20) 



yr{T*) - y,{T*) = - {{Ma{T*)u^ ,u^) - {Ma{T*)u,,u,)) + — (||<r " Wnof) 



First notice that 



RO II 



ImoII, Va;G[0,27r]. 
Moreover, setting Uq = uq + v'^ we get (we omit T* in the argument of Ma) 

{Mau^,u^) - {Mauo,uo) = {Mav'^^v^) + {{Ma + M'a)v^,uo) 







since the first 2x2 block of AIa that is zero at T* and has nonzero component only in the first 
two entries. 
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Remark 19. Note that (i) and (ii) are just the manifestation of the fact that, forgetting about the 
second vertical component 1/2, we are facing the corank 1 case, for which T* is a cut time and there 
is a rotational symmetry that implies that it is also a conjugate time. 

(iii). Now one can proceed in a similar way and compute 

1 /■* 

y2{t) = - / x{s)'Au{s)ds 



2 _ 

2 

{C{t)uo,uo) , 



f (A~\e'^ - I)yAe'^dsuo 
Jo 

u'o [ {e-'^ -I)A~^Ae'^dsuo 
Jo 



where we set 



C{t) = i / (e"'^^ - I)A-^Ae'^ds. (21) 







Since all matrices appearing in ( pi] ) but A are 2x2 block diagonal, the first 2x2 diagonal block of 

K{s) := {e-'^ -I)A~^Ae'^, 

is the product of the respective blocks. A direct computation shows that it is 

ao A— cos(airs) — sin(airs) \ 
ai y sin(airs) 1 — cos(airs) J ' 

where ( ^ "^"^ | denotes the first 2x2 block of A. Integrating from to T* one obtains for the 
\-ao J 

first block of C(r*) 

As before, we set Uq = uq + v'^ and we get (omiting T* in the argument of C) 

y^{T*) - y2{T*) = {Cv^ ^v^) + {{C + C')v^ ,uo) . (23) 

Using (H) and 

\\v^\\'^ = ((cosw — l)ni + s\iiujU2)^ + (— sinujui + (cosw — l)ti2)^ 
= A{ul + ul)sm^{u}/2), 

one gets that ( p^ is linear with respect to the variables 

cosw — 1 = — 2 sin^(a;/2), sino; = 2cos(cj/2) sin(w/2). 

In other words, if we prescribe that the expression ( |23| ) is zero, we get 

Co sin(w/2)(Ci cos(w/2) + C2 sin(w/2)) = 0, (24) 

for some suitable constants Co, Ci, C2 that do not depend on u. The Claim is proved since equation 
has always a nontrivial solution u) G [0, 2tt]. □ 



12 



Let us now show that 7(t) cannot be optimal after T*. From the previous computation we have 
7(r*) 7^ 7'^(r*). By contradiction if 7 is optimal after time T* then the concatenation of 7'^|[o,t*] 
and 7|[T*,T*+e] (fo'^ some e > 0) is optimal as well, which is impossible since all optimal trajectories 
are projections of the Hamiltonian system associated with (|l^) and they are smooth. 

2.3.2 Optimality of geodesies 

In this section we prove that 7(t) = E^°''^''^{t), r > 0, is optimal up to its first Maxwell time 
T* = 27r/(air). 

To this extent, consider the following auxiliary optimal control problem: 

P. Let T < T* and set {x,yi,y2) = "){T). Find a length-parametrized trajectory of the system 

Xi — Ui, i — 1, . . . , Ul, 

yi = \x'Liu, (25) 
,2)2 = \x'L2U, 

starting from the origin, and reaching the hyperplane {x = x} in time T, maximizing the yi 
coordinate. 

Remark 20. Notice that yi > since r > implies that the trajectory is not contained in the 
hyperplane {yi = 0}. 

Lemma 21. The following assertions hold: (i) There exists a solution 7* of the problem P. (ii) 
7* is a length minimizer. (Hi) 7*(t) = E'^'-'''^'^{t) for some {uo,r). 

Proof. Let us prove (i). In problem P, since we deal with length-parametrized trajectories, we can 
assume that the set of controls in (^) is U = {\\u\\ < 1}. The existence of a solution of P can be 
obtained with standard arguments using the compactness and convexity of the set of admissible 
velocities (see P, p^). 



To prove (ii) assume by contradiction, that there exists a trajectory of (25) reaching the point 
(x, ^1,^2) ~ 7*(^) time Tq < T. By small time controllability there exists a trajectory of system 
(I25I ) reaching in time T the point {x,yi,y2), with yi > yl contradicting the fact that 7* maximize 
the yi coordinate. The fact that 7* is also a length minimizer follows from Remark 

To prove (iii) observe that 7* satisfies the Pontryagin Maximum Principle (see again [^) for the 
problem of minimizing —yi = — i/idt = — x'Liudt, i.e. with the Hamiltonian 

m 

nu = Y^ (A, UiXi) + ux'Liu (26) 

i=l 

= pu + rix'Liu + r2x'L2U + ux'Liu. 

where A = {p, ri,r2) are the dual variables to (x, 2/1,2/2) in T*M. In formula (|2^ ) v \s a. nonnegative 
constant. The Hamiltonian equations give 

Vi = -ff = 0, 

P' = = -inLi + r2L2 + uLi)u. 
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Since the final point is constrained on the set {x = x}, the transversality conditions give ri = 
0,^2 = 0. Hence we have 

'Hu = {p + ux'Li)u, 

p' = -vLiu. (27) 

Notice that actually > 0, otherwise the trajectory is a straight line contained in the plane 



{ui = 2/2 = 0}, see Remark 20. The maximality condition and the condition that the final time is 



fixed in such a way that trajectories are parameterized by length give 

'Hu(t){x{t),yi{t),y2{t),p{t),ri,r2) = maxnv{x{t),yi{t),y2{t),p{t),n,r2) = 1, 

V / y 



Notice that a geodesic for the problem (^) associated with the covector {uo,ri,r2) corresponds to 
a control 

u{t) = p' — riLix — r2L2X, (29) 

where 

p' = {-riLi - r2L2){p' - riLix - r2L2x). (30) 

Comparing equations ( |27| ) - (^8|) with (^) - (|30|) it follows that 7* is a geodesic for the problem 
( p5| ) corresponding to an initial covector (uo,z^, 0), for some uq, with ||So|| = 1. Then (iii) is proved 
for r = u. □ 

We have the following 

Claim. 7* = 7, i.e. uq = uq and r = r. 

Proof of the Claim. It is enough to prove that the parameters U(),r such that a geodesic 7(t) = 
{x{t),yi{t),y2{t)) = E'^°'^'^ satisfies x{T) = x with T < T* are unique. 
From the computations in Sections we know that 

x{t) = MAit)uo, where A = rLi, and MA(t) = A-^^e^"^ - I). 

In particular, using the non singularity of A, the equality at t = T gives 

uo = MxHT)x. (31) 

Computing the norm of vectors in equality ([3l|), it follows 

l = IKf = E 'if 4;Jt/2) ' "^^^^ P.(x) = (4-i+4)^/^. (32) 



Notice that the right hand side of ( |32D is the sum of monotonic functions with respect to the 
variable rT, on the segment [0,27r/ai] (T < T* implies rT < 27r/ai). 

Moreover since the curve is length-parametrized we have ||x(T)|| < T. As a consequence there 
exists a unique solution rT of equation (32) in the segment [0,27r/ai]. In particular r is uniquely 



determined and uq is uniquely recovered from equation (pl|). □ 
Since 7 = 7* and 7* is length-minimizer for every T < T*, it follows that tcut = T* . 
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2.4 First conjugate time 

In this section we prove that in the corank 2 case, the cut time is not equal to the first conjugate 
time, in general. This is deeply different from the corank 1 case, where the cut locus always 
coincides with the first conjugate locus. 

It is enough to show that the cut time is not conjugate in the (4, 6) case. Define the Jacobian 
of the exponential map 



Je{t,Uo,ri,r2) := det 



dS d£ d£ as 
dt ' duQ ' dri ' dr2 



(33) 



Remark 22. Recall that the first conjugate time tconj for the geodesic corresponding to the covector 
(zio, ri,r2) is the first time t > for which we have 



j£(t,Uo,ri,r2) = 0. 



(34) 



We have to prove that equation (|34D is not satisfied when t = tcut- 

To compute J^- we use the following trick. Let £ = pdx + ridyi + r2dy2 be the Liouville form. 



Lemma 23. We have £ — - 

V dt 



1, £ 



d£ 
duo 



Proof. The first equality follows from the fact that the Hamiltonian is homogeneous of degree 2. 
Indeed set A = (p, ri, r2) and q = {x,yi,y2), we have £ = Xdq and 



-1) 



A, 



dq 

m 



1, 



since length-parametrized trajectory belong to the set {H = 1/2}. The second and the third 
identities follow from the fact that the Liouville form is preserved by the Hamiltonian flow, hence 
the values of £(^^) and £(^) are constant with respect to t. In particular at t = they are 
annihilated by the Liouville form. □ 

If we compute the exponential map in a neighborhood of a geodesic with (ri,r2) = (r, 0), with 



r 7^ 0, using the identity ridyi = £ — pdx — r2dx2 and Lemma 23 we get 



, fd£ d£ d£ d£\ 
Js = dxA dyi A dy2 ^ 
\ dt duo dri dr2 ) 

1 , ( d£ d£ d£\ 

= —dx A dy2 ^ , 

ri \ ouo dri or2 I 



(35) 



where £{t,uo,ri,r2) = {x{t, uq, ri, r2) , y2it, uq, ri, r2)) denote the exponential map where yi is re- 
moved. More precisely, (p^) is the function of {t,uo,ri,r2) given by 



■det 



/ dx dx dx dx dx \ 
duo duo duo dri or2 



n 



dy2 



\duo ^ 



dy2 dy2 
duo duo 



dm dm 
dri dr2 / 
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where vi,V2,V3 are 3 independent tangent vectors to the 3-sphere {uq G R^, ||uo|| = 1}. We 



Ui 



V / 



, V2 



( \ 



v^3 y 








Prom the computation of Section 2^ one easily gets 

dx 



dun 



dx 
dri 



dx 
dr2 



Moreover (see again Section 2.3) 



where 



■ = MA(t) = A-^(e*^-I), 

-A-^A{A'\e'^ - I) +te'^)uo 
-iMA{t) + te'^)uo, 

-A-^AiA-^e'"^ - I) +te'^)uo 
-A-^A{MA{t) + te^^)uo. 

y2{t) = {C{t)uo,uo) , 
1 



C{t) = — / (e"'^^ - I)A^'Ae'''ds. 
2 Jo 



The function y2{t) from (|37|), being a quadratic form with respect to uo, gives 

'-Vi = {{C{t) + C'{t))uo,v,). 



dm 

duo 



Now we compute these derivatives at i = tcut = ^ where a > b are the moduU of the eig 
of ^. 

It is easily seen that 



B := MAitcut) 



fO 




\0 








\ 





sin(27rfe/a) 2 sin^ (vrfe/a) 



br 



br 



2sin^(7rb/a) sin(27r6/a) 



br 



br 



from which it follows that 



[Bvi Bv2 Bv^ 












n 2sin^(7rfe/a) 
^3 



\0 Ui 



2sin^(7rb/a) 
br 



uA 



+ ^^3 






sin(27r6/a) 2 sin-^ (vrb/a) 



br 



sin(27rfe/a) 
br 



Ui- 



Ui 



br 



sin(27r6/a) 
br 




Mj ' 
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where the last identity defines the matrix M. 

The Jacobian determinant of the exponential map computed at t = tcut is then expressed as 
follows 



Js{t 



cut J 



det 



det 



det 



/ dx 
^ — "1 


dx 
duo 


dx 
duo 


dx 
dri 


dm 

\duo 


dy2 
duo 


dy2 
duo 


dm 

dri 


/ : 










BV2 




dx 
dri 


dy2 
\duo 


dy2 

^^2 

duo 


dy2 
duo 


dy2 
dri 


( ° 
























Mil 


M12 





M21 


M22 



dx \ 
dr2 

dm 
dr2/ 

: \ 

dx 
dr2 



\{Cuo,Vi) {Cuo,V2) {Cuo,V3) ^ 



dxi 
dri 


dxi \ 
dr2 


8x2 
dri 


dx2 
dr2 


dx-j, 
dri 


dx3 
dr2 


dx4 
dri 


dx4 
dr2 


dy2 
dri 


dy2 , 
dr-i/ 



where we use the notation 



M 



Mil M12 
M21 M22 



From (p9D it follows 



where is the matrix 



Js{t 



cut I 



{Cuo,vi) • detM • detiV, 



9x1 


dxi 


dri 


dr2 


dx2 


dx2 


dri 


dr2 



(39) 



(40) 



It is easy to see from the explicit expression of the geodesies that, in the general case when 
a 7^ 6, the three factors in (^) do not vanish identically in no, since the matrix A is arbitrary. The 
proof of Theorem ^ is then completed. 

Remark 24. Notice that M is the zero matrix when a = b. Hence, in the (4, 6) case, t^ut = tconj for 
those 9 such that Lg has double eigenvalue. Moreover in this case the rank of the Jacobian matrix 
drops by 2, since the first three columns are proportional. 
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3 The nilpotent (4, 6) case 



In this section we restrict to the (4,6) case. By the previous discussion the geodesies of the sub- 
Riemannian metric can be written as follows 



' x{t) = A-\e'^ - I)uo 



yiit) 



2./0 



x{syA u{s)ds, 



y2{t) = I x[s)' Au{s)ds 



and we can assume the matrix Lq to be diagonal 



while Lq is an arbitrary skew-symmetric matrix 



Le 



"0 





A = rLg, A = rLg 



/ a 0\ 

-a 

6 

V 0-6 0/ 



/ 

-ao 

-ai —as 

\—a2 -Q!4 







(41) 



a >b, 



(42) 



as 



-as 0/ 



3.1 Proof of Theorem 11 

Recall that the cut time tcut coincide with tcon j if and only if tcut is a time that satisfies the equation 



j£{t,U,9)\ ^ = 

{P2) => {PI). We consider separately the two cases: 



0. 



(43) 



(a) 



Li,L2 both belong to the same subspace, either Q or Q. Then it is not restrictive to assume 
that Li,L2 G Q. In this case all linear combination of Li,L2 belong to Q, i.e. Lq = 
cos(0)Li + sin(^)L2 € Q for every 9. In particular Lg has a double eigenvalue for every 
6 € [0,27r]. Prom the computation of Section § it is easily seen that a = b implies M = 0, 
hence from (^) it follows that tcut = tr 



(b) Li G Q and L2 & Q {Li and L2 plays the same role). By (55) we have [Li,L2] = 0. Let us 
prove then that this property implies (PI). 

Indeed every two commuting skew-symmetric matrices can be block diagonalized simulta- 
neously in the same basis. Hence we can assume that, choosing an appropriate coordinate 
system y 1,1/2, that both Li,L2 are diagonal. As a consequence Lg and Lg are also diagonal. 
Moreover from (|3^ ) it is easily seen that, if both Lg and Lg are diagonal, C is 2 x 2 block 
diagonal, with the first block equal to cl, for some constant c (see also (|2^)). 

In particular it follows that {Cuq,vi) = and again ( ^0[ ) implies tcut = tconj- 
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(PI) =^ (-P2). By assumption the identity 



J€{t,uo,9)\^^^^^^^g^ = {C{e)uo,vi)-detM{uo,e)-detN{uo,e) = 0, (44) 

holds for every uo (the horizontal part of the initial covector) and every 6. Since the exponential 
map is linear with respect to uq in the x-variable, and quadratic with respect to uq in the Ui- 
variables, it follows that (|4^ ) is an analytic expression of {uo,9) (it is polynomial with respect to 
Uq and trigonometric in 6). In particular one of the three factors in (^) must vanish identically. 
Assume that detM(no,^) = 0. Then from the explicit expression it is computed that 

and since a{6) > b{9) by assumption, det M = implies a{9) = h{0), for all 9. 

From this it easily follows that Lq has double eigenvalue for all 9, i.e. if we write 

Le = q{9)+m, 

it follows that one of |lg(0)|| and ||?(0)|| is identically zero (it is not a restriction to assume ||§'(0)|| = 
0). Hence Lg ^ Q for all 9, that implies in particular that Li,L2 € Q. 

It is not restrictive now to assume that a{9) ^ b{9) for some 9 € [0, 27r]. We show that the 
identities 

(a) detiV(no,^) = 0, 

(b) {C{9)uo,vi) = 0, 

both imply that there exists a choice of the coordinates such that Li S Q, L2 € Q. We give 
details only for case (6), the other one is similar. Considering (6) as an equation in the variables 
01,02,03, 04 (the non diagonal entries of the matrix Lq) it is easy to see that the identity (6) can 
be written as an equation 

F{ai,Ui) = 0, 

where F is a quadratic form in the Ui whose coefficients depend linearly on Oj. Since these equation 
should be satisfied for all no = (ni, . . . , ^4), choosing values 

noG {(1,0, 1,0), (1,0, 0,1), (0,1, 1,0), (0,1, 0,1)}, 

one gets the set of 4 linear equations: 

(aa4 + ha\) cos r] + (003 — 602) sin r/ = 0, 

[aaz — ha^) cos -r] — (004 + ha\) sinr/ = 0, 
< (45) 
(aa2 — haz) cos i] + (aa\ + 604) sin r/ = 0, 

{aa\ + 604) cos ri — (002 — haj,) sin 77 = 0, 

where we set 77 = vrft/a, and for simplicity of the notation we denote a = a{9Q), b = 6(^o)- 
It is easy to show, using the fact that a ^ b, that this system has the unique solution 

«! = ... = 0(4 = 0, 
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which means that Li and L2 are both diagonal. Due to this fact they can be written, as pure 



quaternions (see Appendix iA), as a hnear combination of i,i 

Li = ai + ai, L2 = (3i + f3i. (46) 



1 

2/1 \ (a a\ lyi 



Performing the change of variables 

we find a system of coordinates such that Li = i, and L2 = i, i.e. that satisfies (P2). 
3.2 Proof of Theorem [T^ 

In the paper Q it is proved that, if the sub-Riemannian manifold is regular^ with Hausdorff 
dimension f), the Radon-Nikodym derivative of the spherical Hausdorff measure with respect 
to the Popp's measure /i, denoted fsfi, is given by the volume of the unit ball in the nilpotent 
approximation, namely 

fs^.iq) = :r^, Q € M, (47) 

where Jig is the Popp's measure defined on the nilpotent approximation Gg of the structure at the 
point q. Note that fig is the left-invariant measure on Gg that coincides with the Popp's measure 
of the original sub-Riemannian metric at the point q. 

Remark 25. Notice that in our {m,m + 2) case the structure is automatically regular since, by 
assumption, the distribution has constant rank and with one bracket we get all the tangent space. 

Remark 26. In it is proved that fsfi is a continuous function, which is bounded and bounded 
away from zero, in restriction to compact sets. 

Remark 27. For the analysis of the regularity of ( |47| ) it is convenient to parametrize the nilpotent 
unit ball via the exponential map, as a function defined on the whole fiber in the cotangent space. 
In other words we do not restrict to the set {||tioll = 1} and define for every Aq = {uo,r) € 

f(Ao) = vr(e^(Ao)), 

where H is the Hamiltonian defined in ( |l2|) and e*^ denotes the fiow in T*M of the Hamiltonian 
vector field associated with H. Using the homogeneity property H{cX) = c^H{\), Vc > 0, we have 
that 

e-^(sA) = e*^(A), Vs > 0. 

In other words we can recover the geodesic on the manifold with initial covector Aq as the image 
of the ray {tAo, t € [0, 1]} C T*^M that joins the origin to Aq. 

£{tXo) = 7r(e^(tAo)) = vr(e*^(Ao)) = 7{t). 



sub-Riemannian manifold (M, A, g) is said to be regular if, defining the distributions := A, A'^^ 
A' + [A\ A], the dimension of A* does not depend on the point for all i > 1. 
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Due to the previous analysis and thanks to Remark 27 , we can express the volume of the unit 
ball of the nilpotent approximation at a point p as follows 



r-2T[ i-A(e,p) I- 

V{p)= / / / Jeiu,9,r,p)dudrd9, (48) 

Jo Jo Jb 

where is the Jacobian of the exponential map starting from p, expressed in the new variables, 
B = {uq = (til, . . . ,^4), ||mo|| < 1} is the 4-dimensional unit ball and 

27r 

A{6,P) = TTT^- 

max a{Lg[p)) 

The problem of the regularity of the function (|47| ) is then reduced to the regularity of the function 

p ^ V{p), 

where p is a (6-dimensional) parameter. Since the family of sub-Riemannian metrics is smooth with 
respect to p, the exponential map smoothly depends on the parameter p. As a consequence the 
integrand in (|48|), being the Jacobian of the exponential map, is a smooth function of its variables. 

In addition, the function p 1— > A(6,p) is Lipschitz, being the inverse of the "maximum moduli 
of eigenvalues" function, which is Lipschitz (see In particular A{9,p) admits bounded first 

derivative almost everywhere with respect to {9,p). 

Definition 28. Define the following sets 

- S is the set of p such that 3 9 for which Lg{p) has a double eigenvalue, 

- So is the set of p such that 3 a finite number of 9 for which Lg{p) has a double eigenvalue, 

- Soo is the set of p such that V0, Lg{p) has a double eigenvalue. 

Thanks to Lemma |3^, for a generic sub-Riemannian metric, the set of points p G Is a 
union of isolated points. Moreover, due to the expression ( |54|) of the eigenvalues in terms of the 
quaternions given in the Appendix |4.1| , the fact that Lg has a double eigenvalue for all 9 is written 
as ||(/(^)|| = for all 9 (or the same for q). This condition is equivalent to the equation ||5(^)|P = 0, 
that is analytic in 9. In particular this equation, if it is not identically satisfied, has a finite number 
of solution in [0,27r]. 

Remark 29. Notice that the expression (54) for the eigenvalues, provides a crucial obstruction for 
the generalization of the result to m > 4. 

From this it follows that, for a generic sub-Riemannian metric, the set of critical points S is 
the disjoint union S = Sq U Sqo- Moreover the set of points where p A{-,p) is not smooth is 
contained in S. 

Let us write the volume function, depending on the parameter p, as follows 

f-2n rA(e,p) 

V{p) = / / fi9,r,p)drd9, (49) 
Je=o Jr=0 

where we denote by 

f{9,r,p)= / j£{u,9,r,p)du. 
Jb 
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Recall that / is smooth as a function of all its variables, while A{9,p) is Lipschitz with respect to 
the parameters {0,p). In particular it has bounded derivatives. 

We want to prove that V is at any point pQ. To this extent, let us write 

Vip) = / / f{9,r,p)drd9+ / / f{9,r,p)drde. 

Jo Jo Jo Ja{0,po) 

The function 

r27r pA{e,po) 



/ f{9,r,p)drd9, 
Jo Jo 



is always smooth since it is the integral of a smooth function (with respect to p) on a fixed domain. 
Denote now 

/•2-7T rA(e,p) 

W{p):= / f{9,r,p)drd9. (50) 

Je=0 J A{9,po) 

We are left to prove that W is around pq. Notice that, by definition, W{pq) = 0. 

Assume that po ^ ^- Then, since both functions A and / in ( |50D are smooth, W is at pq 
and the derivative at a point p (in a neighborhood of po)i is computed as follows 

-iP)= / ^i9,r,p)drd9+ —{9,p)f{9,A{9,p),p)drd9. (51) 



dpi Jo JA{e,po) dpi Jo dp. 

Assume now that pQ € Sq. The first term in ( ^l|) is continuous. Moreover, since at pq there are 
only a finite number of 9 such that Lq has double eigenvalue we have 

BA 3 A 

j-{9,p) ^{9,po), a.e. 9 G [0,27r]. (52) 

opi p^po opi 

Since ^ is bounded and / is smooth, by Lebesgue's dominated convergence we have that the 
second term is also continuous. 

Finally, consider the case when pQ € Sqo- Since pq is an isolated point, the partial derivatives 
are defined and continuous in A'py \ {po}, where Np^ is a neighborhood of po. We claim that 

dW 

^— (p) — > 0, when p ^ po e T^oo- (53) 

OPi 

Indeed, by definition of Sqo, the cut time A(9,po) coincides with the conjugate time, i.e. it satisfies 
the identity 

Je{u,9,A{9,po),po) = 0, V0e[O,27r]. 

From this it follows that f{9,A{9,p),p) for all 9, that easily implies (]53|). From the fact that 
W is continuous in Np^, in Np^ \ {po} and has partial derivatives tending to zero for p ^ po, it 
follows that W is in Np^. 

4 Appendix 
4.1 Quaternions 

The Lie algebra so(4) of 4 x 4 skew-symmetric matrices is the direct sum 

so(4) = Q © g, 
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where Q is the space of pure quaternions and Q is the set of pure skew quaternions. 
The space Q (resp. Q) is generated by the three matrices z, j, k (respectively i, j, k) 



i = 





-1 





^ 




1 ^ 





-1 


\ 




1 








-1 \ 


1 






















1 


, k = 








-1 














-1 


, j = 


1 














1 








V 





1 


) 




vo 


-1 





) 




V 1 
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-1 









1 





1 






1 ^ 








-1 \ 


1 






















1 










1 
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, j = 


-1 











, k = 





-1 








Vo 





-1 


) 




V 


-1 





) 




V 1 








/ 



If we endow so(4) with the Hilbert-Schmitd scalar product 

(Li,L2) = ■^trace(LiL2), 

then i,j,k,i,j,k is an orthonormal basis. 

The eigenvalues wi, a;2 of A = q + q satisfy: 

-{u;,,2f = {\\q\\±m?- (54) 

As a consequence an element A € so(4) has a double eigenvalue if and only if A ^ Q U Q. 
Also pure quaternions and pure skew quaternions commute: 

[q,q\=0, q£Q, q£Q. (55) 



4.2 Transersality lemma 

Let S be the set of (m, n) smooth sub-Riemannian metrics over M, equipped with the Whitney 
topology. Due to the C°° structure, we have the existence of smooth bump functions and the 
results in this section are essentially local. Then we can assume that S is the set of m-tuples 
F = {Xi, . . . Xjn) of smooth independent vector fields on some open subset M of M", satisfying 

TqM = Aq + [A, A],, for every q G M. 

The vector fields Xi , . . . , X^ form an orthonormal basis for the sub-Riemannian metric g they 
specify. 

Let B be the bundle over M whose fiber at (7 € M is the variety of /c-dimensional vector space 
of g-skew symmetric endomorphisms of Aq. 
Let us consider the mapping 

p: S X M B 

iF,q) ^ 

where Cq has been defined in Remark |l]. It is clear that p is C°° . 
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Let us fix a point Fq £ S, a point qo € M and coordinates (x, y) in M such that qo = (0, 0) and 
the nilpotent approximation of Fq reads in control form 



Ui 



l,...,m, 



yi 



(56) 



In the coordinates y, the space is the vector subspace spanned by the matrices Li, . . . L^. We 
have a natural gradation in formal power series of (x, y) induced by setting that the Xi have weight 
1 and the yi have weight 2. This induces a formal gradation on formal vector fields on Mg^, in which 
^ have weight —1 and ^ have weight —2. The vector fields of the nilpotent approximation ( ^6[ ) 
have weight —1. 

In control form the sub-Riemannian metric Fq itself reads 



/ i \ 

\yk J 



f u \ 

\ \x'LkU ) 



+ H 



(57) 



where H \s a, term of order > — 1 as a formal u-dependent vector field. Then we take a smooth 
bump function b{x,y) which is compactly supported in M and which is 1 in a neighborhood of 
qo = (0,0). We consider the affine space A of variations of Fq of the form 



/ i \ 

m 
\yk J 



/ n \ 

^x'Liu 



\ \x'LkU ) 



+ 



/ \ 

\x'bL\u 



\ \x'bLkU ) 



b{x,y) + H. 



(58) 



This defines new sub-Riemannian metrics Fq + 6F. Since S is open in the set of all rank m smooth 
sub-Riemannian metrics over M, then, for a small perturbation 6F, we have Fq + 5F € S. 

To show that pq^, defined by Pq^iF) := p{qQ,F), is a submersion at qo on the fiber Bq^, it is 
enough to observe that 

Pqo-.A so{4f 

6F ^ {Li + 6Li,...,Lk + 6Lk) 

is an affine submersion. Then we have proven the following Lemma. 

Lemma 30. The map p is a submersion. 

Now let us restrict to the (4,6) case. 

Definition 31. We say that a point go & M is critical for a sub-Riemannian metric F if all elements 
of the subspace £^ (from Remark ||) have a double eigenvalue. 

By formula this means that, whatever the coordinates y, the matrices Li{qo), L2{qo) both 
belong either to Q or to Q. 
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The dimension di of the fiber of the bundle B is the dimension of the Grassmannian G(2, 6) of 
2-subspaces of M^, i.e. di = 8. 

The dimension of the set of pairs Li,L2 that both belong to Q (respectively Q), is the dimension 
d2 of the Grassmannian G(2, 3), i.e d2 = 2. 

Let us define now the partially algebraic "wrong set" W C B as follows: the fiber Wq^ is the 
set of 2-subspaces of the g-skew symmetric endomorphisms of A^j, , whose elements have a double 
eigenvalue. The codimension of in i? is di — d2 = 6. 

The next Lemma follows from Lemma ^ and a non-compact version of Abraham's parametric 
transversality Theorems 

Lemma 32 ((4,6) case). The set of sub-Riemannian metrics that have only isolated critical points 
is open and dense in S. 
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